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Abstract: Leptogenic supersymmetry is a scenario characterized by cascade decays with 
lO copious lepton production. Leptogenic models have striking signatures that can be probed 

by the LHC even in the 10 TeV run with as little as 200 pb _1 of data, provided the squark 
masses are about 1 TeV. Leptogenic supersymmetry spectrum arises in several well-motivated 
models and its signatures are long-lived sleptons, numerous isolated leptons, abundant Higgs 
production, rather energetic jets, and no missing energy. The Higgs can be discovered in 
the h — > bb mode via the 4 leptons+4 jets channel because the leptons accompanying Higgs 
production suppress the background. The superparticle masses in leptogenic supersymme- 
try can be measured efficiently due to lack of missing energy and high lepton multiplicity. 
We estimate that 1 fb _1 of integrated luminosity is sufficient to determine the light Higgs, 
neutralinos, charginos, slepton, sneutrino and squark masses in a 14 TeV run. 



*On leave of absence from York University, Ontario, Canada. 



Contents 



1. Introduction 1 

2. Lepto-SUSY spectrum and decays 3 

2.1 Decay chains 5 

2.1.1 Leptons galore 5 

2.1.2 Higgses galore 5 

3. LHC phenomenology 6 

3.1 Sleptons or muons? 7 

3.2 Production mechanism of lepto-SUSY 9 

3.3 Four-lepton channels 9 

3.4 Five-lepton channels 12 

3.5 Six-lepton channels 15 

3.6 Discovery of the Higgs 16 

3.7 Dirac or Majorana gauginos 18 

4. Models with lepto-SUSY spectrum 19 

5. Conclusion and outlook 21 



1. Introduction 

What makes supersymmetry leptogenic is best explained pictorially in Fig. 1. The decay kine- 
matics in leptogenic supersymmetry (lepto-SUSY) ensures that multiple leptons are produced 
in every decay chain. The production of new particles is dominated by QCD production of 
squarks and gluinos, which are assumed to be at the top of the mass spectrum. Colored par- 
ticles decay into lighter charginos and neutralinos. The charginos and neutralinos are heavier 
than the sleptons and therefore decay into leptons and sleptons. All sleptons decay into the 
lightest slepton which is the next-to lightest supersymmetric particle (NLSP). The NLSP is 
collider stable and eventually decays into the gravitino. Because the gravitino LSP is created 
outside detectors it is not shown in Fig. 1. 

This particular hierarchy of masses is responsible for collider signals with very little 
Standard Model (SM) backgrounds. For example, the leading production mechanisms, qq 
and qq, lead to at least two hard jets from the decays of the two squarks. The two jets are 
accompanied by two sleptons and a handful of leptons as every supersymmetric particle must 
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Figure 1: Lepto-SUSY spectrum and typical decay channels. 



decay into the slepton NLSP. Because of the large mass differences between the squarks and 
sleptons the NLSP sleptons are quite energetic and a large fraction of them have velocities 
larger than 0.95. Some of the fast sleptons would likely be misidentified as muons and we 
refer to the misidentified sleptons as leptons. The total number of observed leptons in an 
event varies depending on the details of the decay chain and how many of the leptons are 
neutrinos and taus, which are difficult to reconstruct. One expects anywhere between two 
and eight observed leptons some of which would be misidentified sleptons, as well as up to two 
properly identified sleptons. We analyze in detail events with 4, 5, and 6 lepton-like particles, 
as events with fewer identified leptons have non-negligible backgrounds and events with more 
leptons are rare. 

The decay chains that pass through the heavier sleptons lead to significant Higgs boson 
production. This happens most of the time for the heavy stau decays and a fraction of the 
smuon decays as long as the mass splitting between the sleptons is larger than the Higgs mass. 
The Higgs is then associated with a clean four-lepton signature and can be discovered in the 
bb channel. This is thrilling as for most other scenarios with a light Higgs, the bb would not 
be the first channel to be observed, and in fact that channel may never be observed. 

Assuming squark masses of about 1 TeV, we show that very significant excesses in every 
(4, 5, and 6) lepton channel can be obtained with only 200 pb _1 of integrated luminosity 
in the 10 TeV run of the LHC. Some of the masses can be reconstructed with this little 
luminosity, while most masses can be determined using 1 fbr 1 of data collected at 14 TeV. 

Our approach is model independent as we parameterize the spectrum and do not make 
assumptions about its origin. However, lepto-SUSY spectrum does arise in several models. 
Low-scale gaugino mediation (LgM) [1] is a class of models with a parametric suppression 
of sfermion masses compared to the gauginos. Other examples include gauge mediation 
(GMSB) [2] with a large number of messengers as well as models with Dirac gaugino masses [3] . 
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The structure of this article is as follows. In the next section, we present a sample spec- 
trum and describe the most important slepton decay channels. In Section 3, we turn to LHC 
phenomenology. We discuss events with multiple leptons, analyze the most promising Higgs 
discovery channel, and study spectrum reconstruction. Section 4 contains our parameteri- 
zation of the soft mass terms. We also outline how the lepto-SUSY spectrum is featured in 
several mediation models. We conclude and outline future directions in Section 5. 

2. Lepto-SUSY spectrum and decays 

The spectrum of lepto-SUSY is as depicted in Fig. 1, namely 

mg,rriq > m^o,m^± > m^ > mh,m^ R . 

The gluino can be either lighter or heavier than the squarks. Such ordering of masses emerges 
in several scenarios of supersymmetry breaking. In Section 4, we propose a natural param- 
eterization of the soft mass terms that is applicable to gauge mediated models. In our pa- 
rameterization, the scalar soft masses are described by four dimensionless numbers and are 
proportional to the gaugino masses. The gaugino masses are assumed to obey the unified re- 
lations so, together with tan/3 and sign /t, we have 7 parameters. The details are not crucial 
for the phenomenological studies we are about to describe and are therefore postponed to 
Section 4. 

We impose constraints on the spectrum that ensure correct electroweak symmetry break- 
ing (EWSB), and that the slepton [4] and Higgs masses are above the direct search limits. 
The Higgs mass is most sensitive to the the magnitude of the stop loop correction. The 
bound rrih > 114 GeV implies a lower bound on rriq, which in the large tan/3 limit trans- 
lates to rriq > 700 GeV. The lightest neutralino Xi can be either Higgsino-like or Bino-like 
depending on the relative sizes of fx and mi, where rrii are the gaugino soft mass terms with 
i = 1, 2, 3 corresponding to the SM gauge groups. If \i S> mi, Xi is mostly Bino-like while it 
is Higgsino-like for ji <C mi. 

Using these results, we find the bounds on the mass ratios that are relevant to the 
production and decay channels we will consider: 

3 > rriq/m2 > 1.1, 5.8 > m^/m\ > 2.2, 
12 > m\jmi > 1, 26 > m 2 /m £ ~ L > 1.9, ///mi > 0.2. (2.1) 

When the constraints are taken into account, the gap between sleptons and gauginos can 
be rather large. The MSSM spectra are calculated using the SUSY-HIT program [5] and a 
sample spectrum is shown in Table 1. The soft masses and our parameters that yield the 
sample spectrum are described in Section 4. 

In our study the trilinear A-terms are set to zero, so the three generations of squarks 
are nearly degenerate and there is little mixing between the left- and right-handed squarks. 
Thus all three generation squarks can be produced at a hadron collider. As a consequence, 
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Table 1: A sample spectrum calculated with SUSY-HIT using input soft terms described in Section 4: 
fi = 294 GeV, A = 0, and tan f3 = 10. All masses are in GeV. 

the Higgisinos with significant couplings to the third-generation squarks only could also be 
produced in cascade decays. Non-zero A-terms contribute the term Av sin (3 to the off-diagonal 
entries in the squark mass matrix, where v = 174 GeV is the EWSB vacuum expectation value. 
Assuming only the stop A-term, At, is large, a sizable mass splitting between the up-squarks 
and stops of order 5m ~ 100 GeV would require Atv/(2mq) ~ 100 GeV in the large tan (5 limit. 
For our sample point, ~ 1 TeV which corresponds to At ~ 1 TeV. This is unnaturally 
large in models that lead to the lepto-SUSY spectrum, hence we neglect the effects of the 
A-terms. 

In the sample spectrum in Table 1, the NLSP is fi. The decay length of the NLSP, 
produced with energy E, in the laboratory frame is (see e.g. [2], [6]) 

\~E 2 /100GeV\ 5 / mn /2 \ 2 

^-^^i^-'v^r) km - (2 - 2) 

where m 3/ / 2 is the gravitino mass. For the typical masses of our model, gravitinos lighter than 
1 GeV are cosmologically safe, i.e. they evade the constraints from overclosure of the universe 
and Big Bang nucleosynthesis, provided that the reheating temperature is lower than about 
10 7 GeV. 

For simplicity, we will assume that m 3 / 2 ~ lOkeV. Then fi has a very long lifetime and 
exits detectors without decaying. The other right-handed sleptons, ejj and will decay to 
fi through In — ► Irtrtr. Due to the small mass difference m e - fi — — m T = 0.6 GeV, 
the lepton and r produced in the decays have insufficient energy to pass pt cuts. Larger 
mass splittings were considered in Ref. [7]. In our case, the decay lengths of e.R and JIr are 
too short to observe an independent track. Even if the decay lengths were long enough to 
observe tracks, the small mass difference means there will not be visible kinks. Therefore, 
the phenomenology will resemble that of the "slepton co-NLSP scenario" of the well-studied 
GMSB point G2b [8, 9, 10]. Note, though, that the ordering of the light neutralinos/charginos 
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and of the heavier sleptons masses is opposite in lepto-SUSY to that of the G2b point. We 
now turn to the decays relevant for collider phenomenology. 

2.1 Decay chains 

In lepto-SUSY scenario, higgsinos and gauginos must decay to the lightest collider-stable 
SUSY particles, £i=e.R, (1%, f%, and hence produce leptons. The number of leptons depends 
on whether the neutralinos go through a short {\ — ► i\) or through a long {\ — > i 2 —> f-i) 
decay chain. Additional leptons and Higgses are produced in the long decay chains. Phase 
space permitted, f 2 undergoes a two-body decay into Z or h. Those channels are closed for 
e-L so the three-body £l — ► li'l'-y decay dominates. For p, 2 t both the two- and three-body 
decays are open. 

2.1.1 Leptons galore 

In lepto-SUSY, the heavier sleptons can only decay into two leptons and the collider-stable 
slepton. Same-sign leptons (SSL) are therefore produced at a higher rate than opposite-sign 
leptons (OSL), which is another remarkable feature of the lepto-SUSY phenomenology. 

In the following, we will present analytical results for the three-body decay widths of 
and JIl. The effect of mixing will be discussed in Sec. 2.1.2 as it is important for the staus 
and smuons. Here we use the notation of Refs. [7] and [11]. 

The number of SSL and OSL is given by 

4 

r SS L,OSL = T{r L - t-Jlli) « ^3 £ Cij I%> 2 \ (2.3) 

where we have adapted the results in Ref. [7). The coefficients are 

cy = S&NjiN&igiNj! + g 2 N j2 ){ gi N^ + g 2 N* 2 ), (2.4) 

where Nij is the neutralino mixing matrix. 

These three-body decays are typically mediated by the Bino and, for Majorana neutrali- 
nos, are suppressed by its mass with Tosl ~ ™~b while T$sl ~ m ~B ■ The Bino mediated 
decays are partially responsible for the SSL excess over OSL in lepto-SUSY. For our sample 
point, where m^ R /m^ ~ 0.44, the SSL final state branching ratio is about 3.6 times that of 
OSL for both selectron and smuon three-body decays. 

2.1.2 Higgses galore 

Even though the decay of the heavier sleptons into the Higgs boson and the light sleptons 
is suppressed by the small Yukawa couplings such decay is competitive with the three-body 
decay channels for the fa, and it completely dominates for the f 2 . Below, we derive the f 2 
two-body decay widths into the Higgs and the Z bosons. Analogous results apply to the 
smuons. The mass eigenstates are defined as: 

(t\ \ _ / cos Of sin 9 f \ ( tr \ 
\t 2 ) \-sin Of cos 8 f J \?l J ' 
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Figure 2: TssL/rosL ratio as a function of /frig . The neutralino mass matrix corresponds to 
our sample point. 



where Of is the mixing angle, < Of < ir, and m?i < rrif 2 . To the leading order, 

. ^m T tan/? 
smt/f « — ^ s- and cost/fWl. (2-5) 

m| — mi 

The two-body decay widths are then 

w~ ~ \ h\ fJ- 2 Vr cos 2 a cos 2 {20 f) 

r(r 2 -> n + /i) = — hyr^rh) , 

\birmf 2 



g 2 sin 2 ( 2 6>f)m 3 
1287T cos 2 0w m2 z 



r(f 2 -> fi + Z) = v ' , r z ) • (2.6) 



where is the weak mixing angle, while the mass ratios are Tf x = m^/m^ and rw^ 
m h/z/ m T2- The dimensionless functions /i,/2 are defined as 



Mr^rn) = ^(1-rl+rl) 2 -r 2 , (2.7) 
/ 2 (r fl ,r z ) = ((l_r 2 i ) 2 -2r|-2r? 1 r| + r|)/ 1 (r fl ,r z ). (2.8) 

In the decoupling limit, where c a ~ sp, one finds that tan (3 drops out of the ratio of the two 
decay widths: 

r(f 2 -> n + g) _ ((1 - rj) 2 - 2r| - 2r|r| + r|) Z^, r z ) 

For the sample spectrum we chose, fi decays into Higgses 53% of the time while for fix, 
including three-body decays, the branching ratio to Higgses is 44%. 

3. LHC phenomenology 

Our analysis of collider signatures was performed using the following software tools. We 
generated the events at the parton level with the Monte Carlo generator MadGraph [12]. We 
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Figure 3: Number of same-sign leptons and leading jet pt distributions for y/s = 14 TeV and 
C= lfb -1 . 

used a modified version of BRIDGE [13] for particle decays to account for the /Ul — JIr mixing. 
Then, we passed the events through PYTHIA [14] to include showering-hadronization effects. 
Finally, we estimated the detector effects following ATLFAST [15] approach. For example, 
to include energy resolution effects we modified the ATLFAST subroutines for the jet and 
electron energy smearing. The analysis of the events is performed with basic parton level r\ 
and pt cuts on jets, leptons, and sleptons. We also imposed jet and lepton isolation cuts. 
We focus on events characterized by 

1. Large cross section production from squark pair production. 

2. Two hard jets (see Fig. 3 for the pt distribution of jets from the squark decays). 

3. At least four lepton-like particles (leptons or stable sleptons). 

We concentrated on these signals because they are practically background free. There are 
no sources of, real or fakes, four-lepton events with at least two SSL and two hard jets with 
the cross sections in the fb range (see discussion in Section 3.3). Our goal is twofold. First, 
we want to show that statistically significant excesses of events can be observed with little 
integrated luminosity. Second, we want to demonstrate that mass reconstruction of several 
states is possible. In order to do that, we examine the events with four, five, or six lepton-like 
particles. Such events probe different decay cascades and are therefore sensitive to different 
intermediate states. 

3.1 Sleptons or muons? 

Long-lived sleptons are a promising feature of some of GMSB benchmark points [8, 16, 17]. In 
lepto-SUSY, the NLSP is a long-lived stau. Mass splittings between selectrons, smuons, and 
staus are rather small, so one cannot observe e# and jlR decays. As we previously mentioned, 
the phenomenology will resemble the so-called "slepton co-NLSP scenario." A novel feature 
of lepto-SUSY is that pairs of sleptons are always accompanied by leptons. We presume the 
existing studies on long-lived sleptons in GMSB should be modified to account for high lepton 
multiplicity in this scenario, but we expect the essential features will remain unchanged. 
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Figure 4: Slepton pt, lepton minus slepton px, and slepton velocity distributions assuming y/s = 14 
TeV and C = 10 fb _1 . 



Heavy, collider-stable particles appear as muons with a delayed arrival at the muon 
chambers. The dominant SM backgrounds are muons from the b and W decays. Cuts on 
the muon pr and isolation requirements greatly reduce these backgrounds. To reduce the 
background from 6-decays, one applies cuts on the tranverse momentum of typically pt >50 
GeV [18]. On the other hand, in our case the slepton is detected in association with leptons. 
This will reduce b-decay backgrounds and a relaxed pt cut is likely to suffice. Top decays 
have a higher pr, but at the same time they have a smaller production cross section. The 
left plot in Fig. 4 shows the efficiency is high with a pt > 50 GeV cut. 

For stable particles one can infer their mass by measuring the momentum and velocity. 
The momentum is measured with the tracker and the muon spectrometer. The velocity can 
be measured using two techniques: time of flight (TOF) from the muon system or using the 
ionization information from the silicon tracker [8, 17] . Cosmic muons are the main background 
for the TOF techique. To reduce this background, TOF measurement is often correlated with 
an independent measurement of ionization in the tracker. At CMS, this approach is suitable 
for velocities 0.6 < (3 < 0.8 [17]. Further refinements could be used to extend the range of /? 
to about 0.9 [19]. At ATLAS, TOF technique may be applicable up to of around 0.95 [8]. 

Sleptons in lepto-SUSY are decay products of heavy squarks. Hence, the slepton /? 
distribution is peaked towards large values. Indeed, most of our sleptons have velocities 
P > 0.95 and would therefore be likely misidentified as muons. To retain as much of the 
signal as possible, we treat the fastest sleptons with (3 > 0.95 as muons while we assume 
that the slower sleptons are properly identified. We also use the muon smearing, isolation, 
efficiencies and charge misidentification parameters in the slepton analysis. Nevertheless, 30% 
(57%) of our sample contains at least one slepton with j3 < 0.8 (< 0.9) and an early slepton 
mass measurement is possible. Only 3% of the events contain 2 sleptons with velocities less 
than 0.8. In the dominant decay channels, as in 2 1 + 2^ + 2 j channel, the total cross section 
of events with at least one slepton with (5 < 0.9 is approximately 150 fb. 

Another challenge associated with slow-moving long-lived particles is the correct assign- 
ment of bunch crossing. The efficiency for assigning the correct bunch crossing for slow 
moving particles decreases steeply with 0. For ~ 0.8 — 1.0, the efficiency is in the 80 — 100% 
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range, whereas for f3 ~ 0.6 it is only 15% [8]. A cut on the slepton f3 > 0.8 eliminates very 
little signal, as illustrated by the right plot in Fig. 4. 



3.2 Production mechanism of lepto-SUSY 

The main production mechanism for lepto-SUSY events is pair production of squarks 



pp -> qq ,qq, qq 

while the pp — > gq channel is suppressed when gluinos are heavier than squarks, see Fig. 5. 
The cross sections for pp — qq and pp — > qq are of the same order for sparticle masses typical 
in lepto-SUSY. In the study we present in the next sections, production cross sections were 
computed within MG at the leading order and they are compatible with PROSPINO with a 
k-factor of 0(1). 




600 800 1000 1200 1400 1600 1800 2000 
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Figure 5: The production cross sections <r(pp — ► qq + qq + qq) (solid blue line) and cr(pp — ->■ qg) 
(dashed purple) at s/s = 14 TcV for mg/mq = 2 using Prospino [20] at NLO. 



3.3 Four-lepton channels 

Events with four lepton-like objects are best suited for reconstruction of masses of the initial 
colored and secondary color-neutral particles — squarks and neutralinos. We refer to electrons, 
muons, and sleptons as lepton-like particles no matter if the sleptons are correctly identified 
or not. The four-lepton events arise from short decay chains in which both neutralinos decay 
directly to heavy stable sleptons Ir and leptons (that is e and fi only) as depicted in Fig. 6. 
Fig. 6 also shows the muon and electron composition of the four-leptons events where muons 
include sleptons. Other diagrams could potentially contaminate this signal. For example, 
production of Xa an d X2 neutralinos decaying into 

xl^Thf R . 

Since Xa 1S wino-like, it does not decay to Ir directly but instead decays to £l and then to 
£r. Those types of events occur at a lower rate than the signal we are interested in, therefore 
they do not pose a problem. 
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Figure 6: Left: Four-lepton channels. Right: Electron (lower-red) and muon (upper-blue) composi- 
tion of the four leptons. 

Events depicted in Fig. 6 do not have any missing energy hence cutting on would not 
affect our signal. Moreover, accurate missing energy calibration will be difficult in the early 
running, especially when hard jets are present in the event. Fortunately, we found that the 
selection of events in terms of the number of leptons is robust under different cuts. As a 
result, we do not impose a $t cu t m our analysis of lepton channels. We will come back to 
the $t cut in the Higgs search discussion in Section 3.6. 

We select four-lepton events using following criteria 

ni = 4 (including sleptons) with |r/| < 2.5, px > 10 GeV 

and parton level isolation cuts ARu > 0.4, AR^j > 0.4, 

rij ^ 2 with \rj\ < 2.5, p T > 15 GeV 

and post-PYTHIA isolation cuts ARjj > 0.4. (3.1) 

With those cuts, the total cross section is 220 (690) fb for ^/s = 10 (14) TeV which corresponds 
to about 45 (140) events in 200 pb" 1 of collected data. Notice very efficient hard cuts on the 
leading jets pt can be applied, see Fig. 3. 

There are sources of SM backgrounds for events with 4 leptons-l- 4 jets. After applying 
the cuts described in Eq. (3.1), we estimated all SM backgrounds are below the fb level. We 
generated the most important backgrounds with ALPGEN [21]: it+jets, W+jets, VFZ+jets 
and QCD jets. Other SM backgrounds like ZZ+jets are tiny, even before we ask for a hard 
cut on the jets- see Fig. 3 in Ref. [21]. 

More relevant backgrounds come from QCD with jets faking leptons. We used a conser- 
vative estimate of 10 -4 probability of that to happen [8]. Although QCD jets have a cross 
section of about 10 8 pb, requiring four fake leptons cuts down their rate to several orders of 
magnitude below 1 fb. Similar fake rate in backgrounds like VF+jets or VFZ+jets leads to a 
cross section significantly below the signal. 

The it+jets background has an initial cross section of about 1 nb. We estimated the 
6-jets produce an isolated lepton in about 5xl0 -3 . For W decaying leptonically and the 6's 
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Figure 7: Left: Neutralino reconstruction from slepton-lepton invariant mass for s/s = 14 TeV and 
C = lib -1 . Right: Dilepton invariant mass distribution. 

producing isolated leptons, we estimate a cross section of about lfb, which can be further 
reduced by applying a cut on the jets pr- 

Depending on whether sleptons are identified or not, we used different strategies to recon- 
struct the neutralinos participating in the diagram in Fig. 6. In the absence of SM background, 
combinatorial background is the main obstacle to the reconstruction. 

Decay products of the squark decays are rather energetic — see Fig. 4 — and tend to cluster. 
Objects coming from a common decay tend to be near each other in the AR parameter space. 
We used this information to reduce combinatorics. For the events with identified sleptons 
((3 < 0.95), we paired a slepton with a nearby lepton through AR^ selection. We followed 
the same procedure to pair the slepton-lepton pair with the nearest jet. For the events with 
sleptons misidentified as muon, we formed dilepton invariant masses and selected dilepton 
pairs with opposite charge and smaller AR«. The AR discrimination is very powerful. Simply 
taking the average of pairing each lepton with any opposite sign lepton would lead to no 
distinguishable features in the invariant mass distribution. 

The slepton-lepton invariant mass distribution fully reconstructs the masses of three 
neutralinos m^o and m^o . We have assumed that enough many slow sleptons are observed to 
establish the slepton mass. The left plot in Fig. 7 exhibits two clear peaks in corresponding 
to Xi an d Xs- The slepton+lepton+jet reconstruction determines m ? -, see Fig. 8. In Table 2 
we present the Gaussian fits of the neutralinos and squark masses. On the other hand, for the 
misidentified sleptons, their energy is taken to be \p\ where p is the three-momentum, which 
is smaller than their true energy E = ^Jp 1 + m 2 . Therefore instead of peaks, the dilepton 
invariant mass distribution shows edge- like structure with end points at about . m - — m~ . 

V x *-R 

Both the shift of the maximum and the smearing of the distribution, apparent on the right 
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Figure 8: Left: Squark reconstruction for y/s = 14TeV and C = lfb . Right: AR distribution for 
the chosen pair (in contour lines) and the rejected pair (in red). 
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Table 2: Gaussian fits for y/s = 14TeV and C = lfb -1 . The last column refers to the number of 
leptons in the channel. For y/s = lOTeV and C = 0.2fb _1 , each channel event rate is reduced by a 
factor of 13. 

plot in Fig. 7, are the result of "missing mass." The slepton mass is not included in the 
calculations of energy whenever sleptons are misidentified. 

Finally, the difference in AR distribution between the chosen dilepton and the wrong 
dilepton combination is apparent in Fig. 8 — correct pairings have lower AR«. 

3.4 Five-lepton channels 

The five-lepton channels, like the ones depicted in Fig. 9, have a long decay chain involving 
X s1 decaying to &l or v. Fig. 9 also shows the muon and electron composition of the five 
leptons. We use this channel mostly for £l mass reconstruction as well as for gaining a less 
accurate estimate of the vl and X s1 masses. 

Events in the five-lepton channel are selected with cuts similar to those used in four- 
lepton analysis and described in Eq. (3.1), except now rig = 5. The total cross section of this 
channel is 137 (426) fb at yfs = 10 (14) TeV and gives 27 (426) events at 0.2 (1) fb _1 . 
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Figure 9: Upper: Five-lepton channels. Lower: Electron (lower-red) and muon (upper-blue) compo- 
sition of the five leptons. 



For reconstruction, we selected slepton-lepton pairs with smaller AIL» and then formed 
slepton+2 lepton clusters by adding to the slepton-lepton pair another nearby lepton chosen 
through AR selection. The decay £l — > 2£ + £r can be fully reconstructed by the slepton+2 
lepton invariant mass. The distribution in Fig. 11 shows a resonance peak at the £l mass, 
fitted to be 252 GeV (see Table 2). 

The v and X s1 cannot be fully reconstructed as their decays involve missing energy. 
Instead, we define the transverse mass variables of the system: 

Mf T = (E lT + $ T ) 2 - (p iT + f T ) 2 , (3.2) 

where i is either slepton-lepton or slepton+two leptons. We cannot estimate $t mis-calibration 
in the early running, and the information obtained from transverse mass, as in Fig. 11, very 
much depends on how accurate the achieved calibration will be. 

The missing energy measured by the calorimeter, Fig. 10, is small and is mostly be- 
low 100 GeV. Conservatively, we impose no missing energy cut since selection on the num- 
ber of leptons is sufficient to reduce backgrounds. Fig. 10 also shows the contours in the 
($ T , Mj; +2iT ) plane. The contours are the most dense around $t ~ 30 GeV, M^ +2£T ~ 
300 GeV corresponding to the Xi ■ Because the missing energy comes solely from the single 
neutrino produced by either v or x^ 1 decay, one expects to see Jacobian edges at M~ i+iT ~ M„ 
and at Mp, 2iT ~ . The plots in Fig. 11 depict an edge at the sneutrino mass 236 GeV in 
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Figure 10: f T distribution and ($ t> Mt) distribution. 




Figure 11: Chargino-sneutrino transverse mass, slepton invariant mass reconstruction for y/s = 
14TeV and C = lib' 1 . 



the slepton-lepton transverse mass distribution and two edges at the chargino masses 294 GeV 
and 677 GeV in the slepton+2 lepton transverse mass distribution. 
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3.5 Six-lepton channels 

Channels with six leptons and no missing energy arise mostly from the combination of two 
types of cascades 

Short: x° -»■ £~r£, (3.3) 
Long: x° -» h I -> £'r £' 2£ . (3.4) 

In terms of the branching ratios, the dominant cascades are combinations of the long (L) and 
short (S) cascades. The (S,L) cascade combination originates from the following gauginos 

^^),(XS^X2).(^X8),(^,XS), (3-5) 

which are ordered by the total production rate from the highest to lowest. The first neutralino 
in any bracket yields a short cascade, while the second neutralino a long one. 



L = 1 fb 1 




Figure 12: Left: Six-lepton channel. Right: Electron (lower-red) and muon (upper-blue) composition 
of the six leptons. 

There are also some (L,L) combinations producing six leptons and no missing energy, 
where a tau is identified as a jet. Among those, only one combination {xa, Xi) 1S produced at 
a comparable rate to that of (x2jX;j)j which itself is very rare. Therefore, our signal mostly 
comes from a combination of the (S,L) cascades. 

Events in the six-lepton channel are selected with cuts similar to the ones described in 
Eq. (3.1), except now ri£ = 6. The muon and electron composition of the six leptons is shown 
in Fig. 12. The total cross section of this channel is 70 (225) fb at ^/s = 10 (14) TeV and 
gives 14 (225) events at 0.2 (1) fb -1 . The slepton-lepton invariant mass would reconstruct 
the neutralino of the short chain with a contamination from partially reconstructing the £l of 
the long chain. See Fig. 13. Those neutralinos were already reconstructed in the four-lepton 
case with better statistics in that channel. 

The slepton+2 lepton invariant mass would reconstruct the long chain three-body decay 
originating from £l. The peak corresponding to the slepton mass is clearly visible in Fig. 13. 
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Figure 13: Six-lcpton channel invariant mass distributions for ^/s = 14 TeV and C = 1 fb 

Using again the AR selection, one can pair the combination of the slepton+2 lepton and 
nearby lepton from the long cascade. The long cascades reconstruct mostly the X3 an d x\ 
since these neutralinos are the most likely ones to originate the long cascades, as described 
in Eq (3.5). Fig. 13 shows clearly the peaks for the X34 neutralinos. In Table 2 we quote the 
Gaussian fit for the Xi mass. 

Finally, we can pair the neutralino reconstruction with the nearby (in AR) leading jet. 
One can use either the long or the short combination. The squark reconstruction is similar 
to the four-lepton case but with lower statistics, hence we do not present it here. 

3.6 Discovery of the Higgs 

As in any MSSM model, the SM-like Higgs boson in lepto-SUSY is too light to decay into 
two W bosons. The Higgs then decays predominantly into bb with a branching ratio of about 
80%. Due to large backgrounds, the Higgs searches at low mass are focused on a rarer but 
cleaner decay, h — > 77, which greatly limits their statistics. 

In lepto-SUSY, in contrast to common-lore Higgs searches, the Higgs discovery channel 
is h — ► bb. Because the Higgs is produced in cascade decays it is free of SM backgrounds (see 
discussion in Section 3.3 for details on SM backgrounds). It is therefore possible to discover 
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Figure 14: Relevant process for Higgs mass reconstruction. 



h through the analysis of a clean bb invariant mass distribution. The relevant decay chain is 
illustrated in Fig. 14. 

Energetic jets, including 6-jets, are produced in the final state in several ways, not only 
through the decay of the Higgs boson. For example, fa or T2 at the end of the chain can decay 
to the Z, which in turn can decay to bb. For our sample point, we found that the decays to 
h are competitive with the decays to Z, with branching ratios: 

BR(/2i h(Z) + fa) = 44.1% (35.1%) (3.6) 
BR(fi -» h{Z) + f 2 ) = 53.3% (46.6%) . (3.7) 

The Higgs and Z are decayed inside PYTHIA [14]. 

To reconstruct the Higgs mass, we selected events characterized by 

n\ < 4 with \rj\ < 2.5, pr > 10 GeV 

and parton level isolation cuts ARn > 0.4, AR^j > 0.4, 
rij ^ 4 with \rj\ < 2.5, p T > 15 GeV 

and post-PYTHIA isolation cuts ARjj > 0.4. (3.8) 

We ordered the jets according to decreasing pp. With no b-tagging, we assume the third and 
fourth jets come from the end of the decay chain and therefore are the b-jets we are interested 
in. We then construct dijet invariant mass of the third and fourth leading jets to reconstruct 
the Z and Higgs masses. In Fig. 15, the Z and Higgs peaks are clearly visible with bin sizes 
either of 5 or 10 GeV. In Fig. 16, we show how relaxing the $x cut slowly introduces more 
combinatorics but still retains proportionately the Z and h signals. We chose a cut $t < 40 
GeV, which we consider a conservative choice for early running. The total cross section for 
those events is 100 (320) fb at 10 (14) TeV, which leads to about 20 events in the 200 pb _1 
run at 10 TeV. 

We can use rough cuts on the invariant mass to estimate the significance of these events. 
With 1 fb -1 of data, there are 50 events for mdijet between 60 and 100 GeV and 37 events 
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Figure 15: Invariant dijet mass distribution with $ T < 40 GeV without b-tagging for y/s — 14TeV 
and C = lib -1 . 

between 100 GeV and 130 GeV. But that counting does not take into account two important 
facts. First, the Z and h peaks overlap. We used Gaussian fits to estimate that the Higgs 
peak contains 51 events and the Z contamination under the h Gaussian is 16 events, leading 
to 35 Higgs events at 1 fb _1 (see Table 3 for details). The second important effect is the 
combinatorial background which lies below the two peaks, and extended in a larger range 
of dijet masses. The combinatorial background increases by relaxing the $t cut and could 
be fitted with a parabola or a gaussian. Our conclusion is that a proper treatment of the 
combinatorial background would require full simulation. 





mass (GeV) 


width (GeV) 


Events 


z 


82 


23 


82 


h 


113 


15 


51 



Table 3: Fits of the Z and h peaks with $ T < 40 GeV, for y/s = 14TcV and C = 1 fb _1 . 
3.7 Dirac or Majorana gauginos 

Lepto-SUSY spectrum can arise in models that contain an approximate U(1)r symmetry, 
making Dirac gaugino masses necessary. Using the clean lepto-SUSY channels, it is straight- 
forward to distinguish the Majorana or Dirac nature of the neutralinos. For example, a Dirac 
gluino forbids squark pair production channels, like pp — > qq, while pp — ► qq is allowed [3, 22]. 
In the cleanest four-lepton final state in lepto-SUSY we have simulated, the selected events 
have both same-sign and opposite-sign leptons. By simple combinatorics, both configurations 
have the same probabilities. In the Dirac gaugino case, only the opposite-sign lepton config- 
uration remains. The lepton charge distribution gives a straightforward discriminant of the 
two types of gaugino soft mass. 
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4. Models with lepto-SUSY spectrum 

In this section we propose a model-independent parametrization of the lepto-SUSY soft mass 
terms and discuss models that lead to lepto-SUSY spectra. Our parameterization is general, 
but it is most useful for models in which the supersymmetry mediation mechanism inolves 
gauge interactions. 

We parametrize the soft mass terms assuming that various superpartners acquire soft 
mass contributions proportional to their SM gauge charges. The scalar masses pick up 
contributions from every gauge interaction they participate in, where each gauge group is 
characterized by the corresponding coefficient Ki 

3 

m 2 (R) = Y J C 2 (Ri)K i , (4.1) 
i=i 

where the sum runs over the SM gauge groups and the quadratic Casimir, C2, is (N 2 — 1) /2N 
for an SU(N) fundamental representation while it is 3/5 times the hypercharge squared for 
The R4 denote the representations of the field under the corresponding SM gauge 
groups. Finally, the coefficients Ki's encode the details of the particular supersymmetry 
breaking mechanism in play, 

Ki = —mjnf (i = 1,2,3), (4.2) 

where mi are the gaugino masses. 

With unspecified dimensionless rij's, this parametrization is still completely general sub- 
ject only to the assumption in Eq. (4.1). In specific models, like GMSB [2] with large number 
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of messengers or LgM [1] (see also [23], [24]), it turns out that the parameters n% are 0(1 — 10) 
numbers. Our parameterization stresses the fact that scalar soft masses are typically one- 
loop suppressed compared to the corresponding gaugino masses squared in such models. The 
squares of the parameters rij are smaller than the loop suppression and therefore the scalar 
soft masses are smaller, by a factor of a few, than the gaugino masses. 

One additional parameter needs to be introduced in the Higgs sector. According to 
Eq. (4.1), mjj and ni 2 Hu are both positive. Typically, EWSB is triggered by additional 
contributions to m H induced by the large value of the top quark Yukawa coupling. We 
incorporate such a contribution as a free parameter n^: 

X — 2,2 a 3^t 22 fAO\ 

& = ~ m H d + m Hu = -^m 3 n 4 . (4.3) 

We also assume that the gaugino masses obey the unified relations, that is mi oc ol%. In 
summary, our parameter space is defined by an overall scale, given by the gluino mass, four 
dimensionless parameters which parametrize the details of SUSY breaking mechanism, tan (5, 
and sign ji: 

7713, th (i = 1,2,3,4), tan/3, sign^. (4.4) 

We neglect the trilinear A-terms as they tend to be small in the GMSB and LgM. As we 
discussed earlier, the A-terms would have to be unnaturally large to make a sizable impact 
on the mass spectrum. 

In Section 2 we discussed the lepto-SUSY spectrum and constraints on the spectrum. In 
terms of our parameterization, the mass ranges in Eq. (2.1), correspond to 

Parameter Range 

m [2, 5] 

n 2 [0.5, 6] 

n 3 > 1.8 

n 4 > 1.75 

Therefore, lepto-SUSY spectra are obtained for 0(1) parameters. For example, the sample 
spectrum presented in Table 1 is obtained from the input parameters in Table 4. 



m 3 


2000 GeV 


ni 


4.8 


n 2 


3.9 


™3 


2.2 


rt4 


6.7 


tan (3 


10 


sign fi 


+ 



Table 4: Input parameters that correspond to the sample spectrum in Table 1. 
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The parameters rii are 0(1) numbers in models where there is a suppression of the squark 
and slepton masses compared to the gaugino masses of order of a loop factor. For illustration, 
we briefly examine how such 0(1) values of n% arise in some models. 

In GMSB, the gaugino masses are generated at one- loop level, while the scalar masses 
squared at two loops. Both the gaugino and scalar masses squared are proportional to the 
number of messengers, N m . Parametrically 



that is one obtains lepto-SUSY spectrum when the number of messengers is of order of the 
loop suppression factors. The above equation holds for as well since soft Higgs masses are 
generated at the next loop order compared to the scalar masses. 

In gaugino-mediation [25], the scalar masses are generated radiatively through diagrams 
involving the gauginos. Whether or not the scalar masses are suppressed compared to the 
gauginos depends on the mediation scale. In high-scale models [26], the natural one-loop 
suppression is compensated for by large logarithm of the ratio of the mediation scale to 
the weak scale. In LgM one naturally obtains rij ~ 1. The m# — mass splitting, 
parameterized by 714, is generated by two- loop diagrams involving the stop and is of order 
indicated in Eq. (4.3). 

Another example of models with lepto-SUSY spectra are models with supersoft SUSY 
breaking [3]. In such models, SUSY breaking is caused by a D-term vev in a new U(l) gauge 
sector. Due to unbroken R-symmetry the gaugino masses are necessarily Dirac type. The 
SUSY-breaking D-term couples to the visible sector through higher dimensional operators. 
The lowest dimensional operator responsible for gaugino masses is a dimension five operator, 
while the operator responsible for scalar masses is of dimension ten. The suppression of the 
scalar masses compared to that of the gauginos is therefore natural for such models. Dirac or 
Majorana gaugino masses can be easily distinguished using charge asymmetries as explained 
in Section 3.7. 

5. Conclusion and outlook 

Lepto-SUSY is a well-motivated supersymmetric scenario. The ordering of the spectrum 
ensures, due to kinematics alone, copious lepton production in decay chains. Therefore, 
the signals are very clean and the LHC has a tremendous discovery potential. Another 
interesting feature is that the Higgs boson is produced in SUSY cascade decays, which reduce 
the background so much that the Higgs can be discovered via the bb channel. 

We studied a sample spectrum of lepto-SUSY and estimated the LHC potential of discov- 
ering and reconstructing this type of spectrum. The discovery channels have cross sections 
governed by the QCD production of squarks and are characterized by two hard jets, two 
collider-stable sleptons, and at least two leptons. If the sleptons are misidentified as muons, 
which is likely, the signatures are four or more leptons and two hard jets. A discovery of the 




(4.5) 
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stable slepton is possible in the early running of the LHC. Lepto-SUSY with 1 TeV squarks 
generates sleptons with velocities 0.6 < (3 < 0.8 (0.9) at a rate of 210 (400) fb at 10 TeV, 
resulting in 40 (80) stable sleptons with 200 pb _1 of data. 

The Higgs discovery channel is through the h — > bb decay in association with four leptons 
and two hard jets. We estimate the prospects for Higgs discovery are good with less than 1 
ftr 1 of data at 14 TeV, again assuming 1 TeV squarks. 

Using the four-lepton channel one can determine the squark masses and all neutralinos, 
except the second Higgsino and Wino, in the 10 TeV run with 200 pb" 1 of data. In the five- 
lepton channel — which is the only $ T channel — one could further determine the masses of the 
sneutrino, heavy slepton, and charginos using the transverse mass variables. The six-lepton 
channel has a smaller branching ratio, but it will allow one to determine the Wino mass with 
just 1 fb -1 of data during the 14 TeV run. 

There are several directions worth pursuing in the context of lepto-SUSY. The light 
Higgs is produced in a clean environment, so one may attempt to use it for an extraction 
of the bottom Yukawa coupling. We only studied the light Higgs, but there may be new 
ways to discover the remaining Higgs bosons in lepto-SUSY. We did not investigate how 
changing the spectrum, by altering various mass ratios, affects the signatures of interest. For 
instance, increasing the squark-gluino mass ratio can enhance the squark-gluino associated 
production, which in turn, may offer new signatures and, of course, prospects for the gluino 
mass measurement. We have not fully exploited the flavor information in our analysis. Since 
the collider-stable sleptons will often be misidentified as muons one could take advantage of 
this fact to refine the analysis. Last but not least, for the channels with $x one could use 
more sophisticated kinematic variables to improve the sensitivity. 
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